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Is The Layer C (413-1000 km) Inhomogeneous? 


By 


Haruo Mirr. 


Abuyama Seismological Observatory, Faculty of Science, Kyoto University, 


Abstract 


There are opinions that the part of the earth’s mantle is not homogeneous. gene: 


which was first introduced by BULLEN (1949), 


(BIRCH, 1939) and the density distribution based on WILLIAMSON-ADAMS-BULLEN’s method 
(BIRCH, 1954) have been thought as the evidences to show that these opinions are valid. 
The method developed in this paper shows that the density distribution can be derived 
without the assumption of hydrostatic equilibrium within the earth’s mantle. The results 
obtained in the present paper show that the abnormal characters of the layer C (413-1000 
km) can be explained by the non-hydrostatic equilibrium state of the layer instead of its 


inhomogeneity. 


$1. We calculate the variation of Griinst- 
SEN’S ratio within the earth by considering the 
pressure effect on it alone, as the temperature 
effect on the GrRUNEISEN’s ratio is known to be 
small. The experimental estimation of the 
pressure effect has not yet been performed and 
we can rely only on the theoretical one. 

We can get the relationship between 7¢ and 
V if the relation between y» and V or v and 
V are known, for the definition of the Griinrt- 
SEN’s ratio is as follows: 


=: > 1 
Te Slog V7 (1) 
where 
1 4 
ym=( : N’ ie Z Zs ) (2) 
4n V WY, ae wee 


V: molar volume, 
NN: Loscumipt’s number, 
Up, Vs: velocities of longitudinal and trans- 


1 
Kyoto, Japan. 
dr 
the calculation of velocities of seimic waves 
versal elastic waves, respectively, 
Ym: DEBYE’S maximum frequency, 
ve: GRUNEISEN’S ratio. 
Birc# (1939) has shown that 
V, 3/2 
°= (142 
yoilt fe) (3) 
vy? =vng(1 +2 fit seeped (4) 
m+2 
Vs" =Vs0(14+2f){1+ f (3m+4)} (5) 


where ff: negative strain, 
Uno" 

Me ID, 
Vso” 


Vi, Vro, Vso: Molar volume, velocities of 
longitudinal and transversal waves at 
zero pressure, respectively. 

We can obtain the following expression for 
ye by inserting the above equations into (1) 
and (2). 


il 1 1 
fe=— + 7 > a Ge No 
3 3 3 f } V3 + 2p oy ( LES ) 
Vso V {1+ f(8m+4)} Er A Ey m+2 
1 13m-+ 4 ilmt 10 
aa 2 m-+-2 
| Vso2V {1+ f(3m+4)} Boallbtet) rae 
112+ 10 
1+f 
m+2 


°3m+6-+4f(3m-+4) 
+ oncty/ (14.7 mt) Gar dt Le 
m+ 2 


(6) 


1+ f(372-+4) 
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On the other hand, Bircs (1952) has obtained 
the approximate expression for 7@ as follows. 
Using the expression of ym in terms of adia- 
batic bulk modulus Ks, density and Poisson’s 
ratio o 


eS oe 
Yn= const. ff ) F(o) (7) 
0 


and the following formula for re 


ve 4) log Ym _ (8) 


we can obtain the following equation 


1 1 (OK ° 
=— 9 
Te eae (9) 


provided that the effect of F(o) on re is neg- 
lible and the equality 0 log K;/0 log p=(0K7/ 
OP), is assumed. Rewriting the second term 
of the right-hand-side of (9) by means of the 
theory of finite strain, we get the expression 
for y¢ as the function of f as follows: 


Go ae! 12+49f 
eu ea 
We have obtained the above two relation- 

ships between y7¢ and f but both methods have 

their good points and bad points. Thus, in 
order to obtain 7¢ by means of (6), we must 
know the velocities of elastic waves at zero 
pressure, but we have no method of estimat- 
ing these values without knowing the material 
composing the earth’s interior. On the other 
hand, y« as the function of / can be obtained 
immediately by using (10) but some errors are 
expected in the estimation, because in getting 

(10), we have made several neglects and as- 

sumptions. Both methods have also some 

errors resulted from the fact that the theory 
of finite strain is not so complete at present. 

Furthermore, elastic constants concerning the 

velocities of elastic waves are the adiabatic 

ones but in deriving (4) and (5), we use the 
isothermal ones. 

Fig. 1. shows the relationships between y¢ 
and f by means of (10) and (6) for various 
combinations of vp, and veo. The combinations 
adopted are as follows: 


Té= (10) 


BS 


oof Ol 02 
Fig. 1. Diagrams of GRUNEISEN’s ratio versus 
negative strain. 
Vpo=7.75 km/sec, Vso =4.35 km/sec. 
—-— + Vpo=7.60 km/sec, Vso =4.30 km/sec. 
SES Uno =8.75 km/sec, Vs9=4.77 km/sec. 
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(1) vpo=7.75 km/sec, vso=4.35 km/sec, 
(2) vpro=7.60 km/sec, ¥so=4.30 km/sec, 
(3) Vyo=8.75 km/sec, Vso=4.77 km/sec, 
(4) by means of (10). 

We see from Fig. 1 that the Gritinetsen’s 
ratio derived from (10) and (6) on various 
combinations of vy. and vs. decreases with 
increasing negative strain in nearly the same 
manner. 

$2. In order to obtain the relationships 
between the depths and GriinrisEn’s ratios or 
negative strains, we divide the earth’s mantle 


into many thin spherical shells. We assume 


(1) Seismic velocities at the depth of 33km after BULLEN (1947), 


(2) (3) Seismic velocities adopted by BIRCH (1939) at zero pressure for the region shallower than the 
depth of 600km and deeper than the depth of 700 km, respectively. 
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J =0.000 at the surface of the earth mantle 
which is at the depth of 33km, and put the 
inner boundary of the first spherical shell at 
the depth where / =0.005. We assume in this 
spherical shell that the composition and 
structure of material does not change and 
GrUNEISEN’s hypothesis (SLatErR, 1939) 


__ const. 
acre (11) 
holds, then we can obtain the following 
equation 
1 2 
el *3 as en = vs3 ), 12 
v, a dee (12) 


( aes 
Ups a Vs? i) 


by inserting (2) into (11). On that occasion, 
we put the Drsyr’s maximum frequency vm 
equal to the lattice frequency » (Urren, 1952). 
Suffixes 1 and 2 denote the physical quantities 
at the inner and outer surface of the spherical 
shell. The left-hand-side of (12) can be esti- 
mated as the value of f at the outer and 
inner surfaces of the spherical shell and the 
mean value of y¢ are the known quantities. 
The denominator of the right-hand-side of (12), 
(Vp *+2us~3)., is the value at the outer surface 
of the first shell, i.e., at the depth of 33 km, 
and can be estimated from the observed seis- 
mic velocities. Thus, we can get the value 
of (vp -?+2vs~3); in (12). We can estimate the 
depth of inner surface of the first spherical 
shell provided that we have the relationship 
between the depth and (v»~*+2v;~%) beforehand 
(Fig. 2). Seismic velocities necessary in the 
' calculation are adopted from ButuEn’s book 
(1947). 

The second spherical shell is formed as fol- 
lows. We may consider the strain at the 
outer surface of the second spherical shell is 
equal to the one at the inner surface of the 
first spherical shell, f =0.005, since the strain 
of the material is continuous within the earth’s 
mantle. We put the inner surface of the 
second spherical shell at the depth where f = 
0.015. The value of (vp-?+2v5-*), at the outer 
surface of the second spherical shell is known 
_from Fig. 2. The mean value of ye in this 
shell can be known from Fig. 1 and the value 


30-3 3 
Upt+2Ve (se¢/cm) 


20 


: — Depth (em) 
0 ~ 1000 2000 3000 
Fig. 2. Diagram of v»y-3+2v,~-3 versus cepth. 
f 
2.0 0.21. 
aS 
Ney 
a 


PS 0u0 


aed 


—=» Depth(&m) 
ae ° 1009 2000 3000 
Fig. 3. Diagrams of depth versus Yg¢ and 


negative strain. 
Vpo=t.7> km/sec, Vs9=4.35 km/sec. 
Vpo=/.60 km/sec, Vs 9=4.30 km/Sec. 
—X—X— Vpo=8.75 km/sec, Vso =4.77 km/sec. 
by means of (10). 
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of V;/V2 from (3), for the strain in this shell 
is assumed to be between 0.005 and 0.015. 
Using the values estimated as above, we can 
obtain the depth which corresponds to the 
inner surface of the second spherical shell in 
the same way as in the first spherical shell. 

The boundaries of the spherical shells in 
the other parts of the earth’s mantle can be 
determined in the same manner. This division 
of the earth’s mantle depends on the assump- 
tions that the composition, structure and 
GRUNEISEN’S ratio are constant in each shell, 
but not necessarily the same among them, and 
the strains are continuous at the boundaries 
of the spherical shells. The accuracy of the 
result can be raised if we divide the earth’s 
mantle into narrower strain ranges. Fig. 3 
shows the relations among negative strain, 
GRUNEISEN’s ratio and depth. 

The values of GrtNEISEN’s ratio of some 
minerals are as follows: 


Table 1. GruNEISEN’s ratios of minerals 
Corundum ere Diopside 0.8 
Rutile OG Quartz 0.8 
Periclase 1.6 Anorthite O27 
Pyrite 1.6 Albite O25 
Olivine eZ Beryl 0.4 
Diamond 1.0 Zircon 0.4 


These values are rather small, compared with 
the values near the earth’s surface which are 
shown in Fig. 3. The reason for this discre- 
pancy is not known. 

§3. The ratio of molar volume at the inner 
surface to that at the outer surface of each 
spherical shell can be obtained from Fig. 3 
and (3). (Table 2) The calculations are per- 
formed for the following two cases 

(1) Vs0=7.75 km/sec, vso=4.35 km/sec. 
(4) by means of (10). 

If we assume that no discontinuous density 
change occurs at the boundaries of each 
spherical shell by the change of composition, 
the lattice structure and etc., the density dis- 
tribution, mass and moment of inertia of the 
earth’s mantle can be estimated, taking the 
density at the depth of 33km 933 as a para- 


Table 2. The ratio of molar volume at the 
inner surface to that at the outer surface 


of each spherical shell. 
er a ee 


Depth (km) 
Shell number Vi/V2 ee 
(1) (4) 
| aa 33 
1 0.9852 | te e 
: 0.9710 | 935 | 9296 
: 0.9716 | 355 | 340 
4 0.9721 | EE rice 
5 0.9726 | spo A ee 
6 0.9731 4 Re 
7 0.9736 oe ae 
8 | 0.9740 678 647 
9 0.9745 | 758 726 
10 0.9749 
| 872 797 
11 0.9753 ee Ses 
12 0.9757 1330 «1171 
13 0.9761 ee ae 
$ | 0.9765 | 913 | 1789 
15 | 0.9769 2389 | 2145 
16 0.9772 POE age 
17 0.9775 | Re pare 


meter. The results of calculations are as 
follows: 
(1) Mn=1.1469933 x 10?" gr, 
In=2.0244 033 x 10“ gr.cm?. 
(2) Mn=1.1591\o53-< 104" er; 


Im=2.0473 033 X 104 gr.cm?. 

The mass and the moment of inertia of the 
earth’s core can be obtained, taking the central 
density gc aS a parameter provided that the 
density distribution in the earth’s core can be 
obtained by the WutnuiAmMson and Apams 
method. Thus, the mass and the moment of 
inertia of the earth can be expressed by using 
two parameters 33, 0c. . These parameters can 
be settled by the observed values for the mass 
and moment of inertia of the earth. The 
results of the calculation are as follows: 

(1) 03s=3.46 gr/em? .=12.41 gr/cm3 

(4) 033=3.42 gr/cem? ,=12.44 gr/cm? 
By using the above values and the results in 
Table 2, we get the density distribution within 
the earth as shown in Fig. 4. 
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Density (qm /cm3 ) 


{1.0 


8.0 


5.0 


—=> Depth (fm) 

(0) 2000 6000 

Density distributions within the earth. 

e Vpo=7.75 km/sec, Vs 9=4.35 km/sec. 
by means of (10) 

after BULLEN. 


4000 


Bircw (1954) has shown that the initial 
density at the top of the mantle must be 
greater than 3.7~3.8 gr/cm? in order to obtain 
the reasonable values for the mass and the 
moment of inertia of the core. But his method 
of calculation, 1.e., Wr1L1amson-ApaAms-But- 
LEN’S method, is based on the assumption that 
the earth’s mantle is 

zsothermal, 

homogeneous and 

in the state of hydrostatic equilibrium. 
This high value of the initial density, he 
thought, is one of the evidences that the layer 
C is not homogeneous. The calculation in this 
section, however, shows that the reasonable 
density distribution can be obtained from the 
assumption of homogeneity of the mantle. 

§4, The distributions of the seismic wave 
velocities can be obtained by means of (4), (5) 
and Fig. 3. (Fig. 5) The solid lines in the 
figure show the observed values. The small 
circles after Brrca (1939) are obtained under 
the assumption that the earth’s mantle is 


Velocity (RM/sec) 


IST 


De pth (42m) 
O (000 2000 3000 
Fig. 5. Distributions of seismic velocities. 


@ by the method described in this paper. 
© after BIRCH. 


divided into two layers which have the differ- 
ent seismic wave velocities at zero pressure 
and each layer is 

zsothermal, 

homogeneous and 

in the state of hydrostatic equilibrium. 
We see in the figure that the black points 
calculated by the present method is a fairly 
good approximation to the variations of seis- 
mic velocities with depth. Therefore, the 
variations of seismic velocities in the earth’s 
mantle can be explained by the assumption 
that the earth’s mantle is homogeneous but 
not in the state of hydrostatic equilibrium. 


§5. The pressure variation with depth is 
dP OIP\ GY ff CHEN GhIk 
= — 13) 
dr Cre , dr +(5r), dy \ 


and equals to —gp if the condition of hydro- 
static equilibrium is satisfied within the earth, 
where 
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Fas Me)=an\" oredr, (14) 
7 0 
G: 6.670 x 10-8 dynes.cm?/gr. 
If the temperature increases with depth, 
dP OP\ dV Kr dV 
= — = 15 
eo cy 7 dr Viaas oe 


because (OP/OT);=Crre/V>0 (from GrRUnet- 
sEN’s equation of state) and dT/dr<0, where 
Cy is the specific heat at constant volume per 
mole. The right-hand-side of (15) can be 
transformed into Ky/(7,.—7)-log V./Vi and the 
approximation of Ky by Ks make the numeri- 
cal calculation of this term possible by means 
of Table 2, Fig. 3 and the values of Ks/o 
derived from the seismic data. 


ce 


Depth (Rm) 
3000 


(0) 1000 


2000 


Fig. 6. Distributions of pressure gradient, 
neglecting the temperature effect. 
@ vpo=7.75 km/sec, vs 9 =4.35 km/sec. 
© by means of (10). 


The results of calculations are shown in 


Fig. 6. The solid line in the figure is go 
calculated by means of (14) and Fig. 3. We 
see from this figure that the layer C is not 
in the state of hydrostatic equilibrium. 

§ 6. 
layer, 


where the condition —dP/dv=gop is 


The temperature gradient in the D> 


thought to be satisfied, may be estimated by — 


means of (13). The results of 


calculation | 


show that the temperature gradient in the | 


layer D is about 1.8°/km throughout this 
region. 
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Wave Generations from Line Sources within the Ground. 


By 
Hitoshi Taxpucni 


Geophysical Institute, Faculty of Science, Tokyo University, Tokyo . 


and Naota KosayAsuHt. 


Department of Precision Mechanics, Faculty of Technology, Chad University, Tokyo. 


Abstract 


Wave generations from line sources within the ground are studied. 


Numerical details 


are worked out for the line source of step function type in time. The results to be noted 


are as follows; 


(1) At the “epicenter”, the displacement is a simple pulse followed by a gradual 


decrease to a permanent displacement. 


(2) RAYLEIGH waves appear at a certain distance from the epicenter. 

(3) Permanent displacements are rather large even at points far away from the origin. 
They are comparable with the amplitudes of RAYLEIGH waves at the respective points. 

(4) There is no S phase in seismograms to be obtained for the wave origin of 


dilatational type. 


§1. The propagation of tremors over the 
surface of a semi-infinite elastic solid due to 
a line source at a certain depth below the free 
surface was discussed by H. Lamp (1904), 
H. Naxano (1925) and E. R. Larwoop (1949). 
These authors were interested in the displace- 
ments at points far away from the wave origin, 
and did not study the tremors near the origin. 
-Furthermore, by using the method. of contour 
integral, they got some expressions for the 
displacements corresponding to the P, S and 
RayLEIGH waves. As is understood easily, 
each component wave above cited has its dura- 
tion time. Thus, in order to get the displace- 
ment at any time, we must sum up the 
contribution from each wave. In the course 
of this summing up process, some of the above 
waves may be masked by some other waves. 
In any of the above papers, this summing up 
process has not been made. Thus, in a strict 
sense, we do not know the exact time varia- 
tions of the displacements even at points far 
away from the wave origin. In short, it is 
required in this kind of problem to get simple 
expressions by which the displacements at any 
point and at any time may easily be calculated. 


In the analysis to follow, such expressions will 
be given, and by using these expressions, some 
numerical examples will be worked out in 
detail. 


§2. It was shown by Nakano that for the 
compressional line source 
_O¢ _O¢ 
Mita , eet ’ 


b= %y eM Hy(MR) 
—- eS asians Odh : 


0 


Re=27+(y—f)? ’ 


we eee 
: by, +24 Vy» 


at a depth / below the free surface, we have 


where 


(2.1) 


Wa, y=0, d= 2ibe| BE elpirtee- did, 


is 
OT ie eee ae 
vie, 9=0,)=R |" 7 es Gas. 
where P(E) =(2&?—k?)?—4 a Be’, 
an oS ey Py py en ee ae = i paee 
a (6 h ) B =z k , k Py I Vs 
(ee) 
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itty, aly 


In order to satisfy the radiation condition 
(SoMMERFELD. A: 1949, p, 189) at 2-00, we 
must have 

the real parts of a and f>0. (2.3) 
We can satisfy these conditions by interpret- 
ing the integral expressions for w# and v in 
(2.2) as the limits of the contour integrals 
u(x, y=0, 2) 

_ogpe| CVC—R ipt+ilta-fVEC-h 

2ik | Ce a 


v(x, y=0, 2) 
_ya| 20°? int +i¢w—fVG— 
4 \, FX¢) ‘ a 
F(O)=(20?-—-RBY4C2V C—BV C2 — FB? 

(2.4) 
when semi-circles around six singular points 
+h, +k and +« are made vanishingly small 
and V€?—h? and V€?—k? at the right end of 
ZL are assumed to be positive real. 


e 
* -bA+h+ 46 


In Fig. 2, +* are the poles of the integrands 
in (2.4) for which 
Ficce\=0. (2.5) 
In order to calculate # and v in (2.4) by the 
method of contour integral, we shall integrate 
the integrands in (2.4) on the path Z+(B—F 
»D—>A-E-G-C), where BF and GC are 
parts of a large circle with the center at the 
origin O and D->A- £ is the path on which 
Z=ile—fVOC—h? (2.6) 
in (2.4) is pure imaginary, and on FD or EG, 
the imaginary part of 7¢zx—fVC?—h? is kept 
the same as that at D or &. The points D 
and E are made far away from A in the fol- 
lowing analysis. Since there is no singular 
point inside the path and since the integrands 
vanish on BF and GC, we may replace the 
integrals on the right-hand sides of (2.4) by 
those on G>E—A-D-F. In the following 
analysis, some physical reasons will be given 
why we may omit the contributions from the 
paths GE and D-F. Thus we have 


u(x, v=; t) 


ih /C2— he ipt+i¢e-fVe—We 
= 2ik? CVO— BR iptt+ica-fV¢ 
; ae F(€) dl, 
va, y=0, Z) 

ap 2¢ —k? ipt—i¢a—fVC—h 

a pret Fe) ° ae i ; ; dl. (2.7) 


$3. In order to transform (2.7) into more 
convenient forms for numerical calculations, 
we shall make some mathematical considera- 
tions in the present section. 
(a) At first, it can easily be shown that the 
coordinate € at A and Z in (2.6) there are 
he 


ge 
Se , 


Z=—thr , 
© 


r=Verty, (3.1) 


respectively. 
(b) ¢, V@—f and V@—F on AD ara 
of the forms —(+) +2 (+), where (+)'s are 


some positive real numbers. € on AE is the 


complex conjugate of that on AD, and Ve?— h2 
and V€2—k? on AE are —(complex conjugates) 
of those on AD, respectively. 


(Cc) Z=ie—fVC%—f* on AD and AE can 
be put 
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Z=—ihr cosh ¢ , 0<d<+0, (2) 
where 
h 
C= ; (—2a cosh 6+7f sinh ¢) , 
di =+V C?— hide , (3.3) 
respectively. 
(d) By means of (3.3), (2.7) is transformed 


into 
ula, y=0, p=ae|"iQeow cosh DAG, 
0 


v(a, y=0, p=28|"ivageor-m cosh Hg, 
0 


(3.4) 
where U; and V, are the real and imaginary 


BA Saad 2h Sle cae 


parts o and 
F(¢) 
Soe ae on AD, respectively. 
(e) Putting 
€=pal’, h=pa=" h=pb=" —mpa, 
p Ss 
m= (GRD) 
a 


and assuming S(z) instead of e’”' in (2.1) for 
the time variations of the displacement and 
stress near the origin, we have 


Gn A 9 Gar cosh O06. 
0 


p=27 | Vue Star cosh ¢)d¢, (3.6) 
0 

‘where S’ means the differentiation with respect 
to the argument in the bracket. 

(f) The radial displacement wr and radial 
stress Ter in the immediate neighbourhood of 
the wave origin are given by (2.1) as 
_S(Z) 

R : 
In the following numerical calculations, we 
shall assume 


Kadai ) 


UR 


,  Trr= SU), 


SH=0,.—o<f<.0 , S()=const.,. “7=—0 
or 
S(H=0, Z<0, 
S(j=a0, S(t dt=S, say at t=0. (3.8) 


This is nothing but the time variation assumed 


in Lapwoop’s paper. 
(g) For the time variation in (3.8), we have 


ule, y=0, Dade =I Ui(’) ) 


ig Fi sinh ow) cosh b= = ie 
v(x, y=0, t)=2m?> Asie, 
Ji A sinh b cosh b= : 
Of r/V ’ 
0<¢<-~, (3.9) 
where 
b rd 
m= en a , v= V 224+ f? ; 


¢=" (—wcosh +if sinh d) (3.10) 
and UVC’) and VC") are 


the real part of A A Serle 
me) 
—1 (26’?—m?) 
ree 


LF 
the imaginary part of C! 


(3.11) 
on AD, respectively. In (3.11), F(€’) denotes 


F(€’)=(26"—m?y?—467V €?—-1 VE 
(Gol) 


and it is to be noted that the real and imagi- 
nary parts-of VC?—1 or Vf?—m' on AD 
are negative and positive, respectively. (3.9)- 
(3.12) are the formulas by which we can cal- 
culate the horizontal and vertical displacements 
on the free surface for the wave origin (3.7) 
and (3.8). The displacements are calculated 
for the non-dimensional independent variables 
x t t 

FMS ae IVs 
distance referred to the depth of the wave 
origin and the time referred to the time of 
transit of P wave across the distance 7= 


which are the horizontal 


eee: 2 
Vz +f : 
§ 4. In our numerical examples, we assumed 
i a SA NI ae: (4.1) 
Vs 


and calculated the displacements at =0, 1y 


a 


De Bie Tah MO), PAO) eh aral 10: if wu and z v thus cal- 


Ww Te) 
culated are shown in Fig. 3. 
The results to be noted are as follows: 
(a) As it should be, there is no displace- 
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a) «/f=0 b) «/f=1.0 
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f) x/f=10 


g) xif=20 h) x/f =50 


Fig. 3. thich line: 


—$—$—_————— eT 


Wave Generations from Line Sources within the Ground. 11 


ment before the time ¢= 7 ~ 


vp 
(b) Horizontal displacement « at x=0 is 
zero all through the time. This will be under- 
stood by the consideration of symmetry as 


shown in Fig. 4. Vertical displacement at 


| . 


GW 
SSS 
| 
1 
1 
vy? 
a | 
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x=0 is a sudden upward motion corresponding 
to the P phase followed by a gradual decrease 
to a permanent displacement. 

(c) Vertical displacement v at a=/ is very 
similar to that at v=0. wz at this point is a 
sudden outward motion corresponding to the 
P phase followed by a gradual decrease to a 
permanent displacement. There is neither S 
phase nor Ray.ercH phase in the displacements 
at =0 and /. 

-(d) In the displacements at «=2/, there 
appears a phase corresponding to the RAYLEIGH 
The time in which the RAYLEIGH wave 


wave. 
5 5 wv v 
traverses the distance x is f= Ve = 1.88 V,* 
1.887, . Corresponding to this time there is 
p 


a maximum in the vertical displacement v. In 
view of these results, we may say that the 
RAYLEIGH wave appears at about x=2/. 

(e) There is no S phase in any of our 
figures. This seems to be inconsistent with 
the result obtained by Larwoop. Thus, in 
Fig. 22-(a) in Lapwoon’s paper which concerns 
with the same problem as ours, we see a 
phase corresponding to the S wave. This con- 
tradiction, however, may be settled as follows. 
As was said in section 1, Lapwoop obtained 


the expressions for the component displace- 
ments corresponding to the P, S and RayLEiar 
waves. It is these component waves which 
are shown in his figure above cited. In order 
to get the total displacement at any time, he 
must sum up the contribution from each com- 
ponent wave. He did not make this summing- 
up process. On the other hand, as our method 
of treating the problem is entirely different 
from Lapwoop’s, it is not necessary in our 
analysis to make the summations. Thus, it 
may be that Lapwoon’s S phase is masked by 
his RAYLEIGH wave, say, in the course of the 
summing-up process, and on making the sum- 
mation, Larwoop’s results agree completely 
with ours. In order to ascertain this point, 
using his formulas for the component waves, 
we calculated the amplitudes of his S and 
RAYLEIGH waves at the time when the S com- 
ponent wave takes its maximum amplitude. 
The amplitude of the RaynercH wave at this 
time turned out to be about 10 times that of 
the S wave. Thus our inference was shown 
to be correct. 

(f) Putting 7, 2z>co in (3.9)-(8.12), we get 
the following expressions for the permanent 
displacements at points far away from the 
wave origin 


a 


S) 


v0. (4.2) 


tee 
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These asymptotic values are shown in the 


figures for 710, 20 and 50 in Fig. 3. From 
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these figures, we see that around these 
asymptotic values, the RayLEicH waves take 
their typical forms as shown in Fig. 5. As to 
the asymptotic values themselves, they must 
coincide with the statical deformations due to 
a statical compressional line source at the 
depth f. This kind of problem was studied 
by B. Jerrrey (1921) using bipolar coordinates, 
and our results can be shown to agree with 
his results. Anyway, it is rather surprising 
that the permanent displacements are of the 
same orders of magnitude as the amplitudes 
of the RayuericH waves for the step-function 
type time variation at the wave origin. 

(g) wand v in (3.9) at a certain ¢ denote 
the elementary waves which arrive at the 
observing point at this time. Since ¢ in (3.9) 
is between t=) and t=co, (3.9) covers the 
whole plausible time interval. Furthermore, 
as was said before, the permanent displace- 
ments in (4.2) were shown to agree with the 
statical deformations for the statical compres- 
sional origin. In view of these results, it may 
be justified to omit the contributions from the 
paths GF and D-—F in getting (2.7). 

(h) Recently, using agar-agar as an elastic 
medium, K. Kasawara (1955) has made the 
experimental study of the present problem. 
His results agree faily well with those obtained 
above. 

§5. In the similar way as above, we can 
study the surface displacements due to a 
distortional line source. It was shown by 
H. Nakano (1925) that for the distortional 
line source 


Og Na 0d 


aia ’ Vo= Ox ’ 
yam en eR) = Aes eonh 6) ad 
0 
where 
Ract+y-f, k=p/ = 6.1) 
pp Vs 
at a depth f, we have 


207 —R? pipt+ia—-fVO~ ke 


- are 
Pee ace) i 


ne, y=0,0)=k| 


v(a, y=0, #) ay 
aie EVE ies tte FVM, 
\, FO) i) 

F(Q)=Q202-—RY-4C2V C2— WV C2— Fk? , (5.2) 
where Z denotes the path shown in Fig. 2. 
Deforming the path of integration as in 
section 3, we get the following results which 
correspond to (3.7)-(3.12). 


2uS(t 
uf R—0)=— » Ter(R-0) a (5.3) 
S®=0,. —0o<#<0; S()=const., #20. 
(5.4) 
S@)=05 205 
Shi=e;, SbHdt=S', 7=0. (5.5) 
ie Jy t yaa 
sei Wa a te = Ca hee. c 
br The z 
wa, y=0,t)=—4 Sf | 
fy \sith¢7 3, 
br "Irs 2 
0<d<~, (5.6) 
¢’=1(—xcosh¢+if sinhg), r=V224f?. 
(5.7) 


Vi(€’) and U, (€’) in (5.6) are 
the imaginary part of Catt AES 


BCC) 
laf 72 LL pan 
the real part of © Y¥S?—-1VE?—m” 
r FC) 
(5.8) 
on AD, respectively. In (5.8) 

F(0’)= (267-1? -467V €2—1 VE2—m? , 

m = 1 eee Us (5.9) 


m b Vo 
and it is to be noted that the real and imagi- 
nary parts of V€?—1 or V€"?—m? on AD 
are negative and positive, respectively. In the 
case of the distortional origin, there is one 
more point to be considered. When 


kx ' : 3 134 moe Vsf 

TVS Ve Veta 
the point A comes to the left-hand side of € 
=—h as in Fig. 6, and the singular point ¢ 
=—h is included in the path of integration 
L>B>F>D>A->E>G>=C. In order to get 


rid of the difficulties thus caused, we shall 


(5.10) 


Wave Generations from 


late ie 


consider the integral along the path L>B—fF 
—D-—A,—H->A,-E-G-C in Fig. 7. Since 
there is no singular point inside this new 
path, we can rewrite (5.2) as follows. 


uevy=0,0=R(| +| -| ) 
Haedyo. ) Ale Ay 


26° —h? pt + ila-fVO— hye 
Ff) 


v(e, y=0, t= —2i84( | +| _| ) 
HA, +AQD A,H AyH 


CV CAs Hi ipt +i¢e-fVeE-kh We: 
rae) 
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The integrations along the path LA,+A,D 
can be reduced to the same formulas as in 
(5.6)-(5.9) and those along the path A\W—A,H 
are transformed into 


u(x, y=0, t)(additional) 


— Bp SL BOP mE) Eye 
sou G¢é) : 
v(z, y=0, t)(additional) 
wt oe — 1) "Ce mp 
— =4) 2 
Re G() 
G(E)=[(2E°— m?)* + 16E*(E?—1)(m?—€?)] 
i al 

iy ks pas __ £2)1/2 

(E21)? , 2 &2)1/2— 4+ real, (5.12) 


The time at which these additional displace- 
ments begin are given by 
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~ Vs cos 9 Vow oon Vile Vee 


sino = | 

Dp 

This is the time corresponding to a kind of 
surface P wave which occurs at point A in 
Fig. 8. From Fig. 8, it can easily be under- 


(5.13) 
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stood that we must have the condition (5.10) 
for the occurrence of this wave. 
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§6. Using (5.6)-(5.9) and (5.12) and putting 


1 
3 


as before, we calculated the displacements at 
i =0,. 1-2, 3, 4.5 and 10. 
obtained are shown in Fig. 9. The results 
to be noted are as follows: 

(a) Vertical displacement v at x=0 is zero 
all through the time. This will be understood 
by the consideration of symmetry as shown 
in Fig. 10. Horizontal displacement u at 2=0 


A=, n@=3, me= (6.1) 


The results 
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is a sudden inward motion corresponding to 
the S phase followed by a gradual decrease 
to a permanent displacement. 

(po) At 2w=f, we can see the surface P 


wave corresponding to (5.12) and (5.13). The 
time Z = E for this wave is smaller than 
br r/Vs 


1. There is no RayzeicH phase in the dis- 
placements at this point. 

(c) In the displacements at «>3/, we can 
see the Ray.ercH phases. The time in which 
the RayLreicH wave traverses the distance x is 


He G y : : 
=~ =],09~-=1.09—. Corresponding to this 
mie 7, Y, PORE, 
time - : =1.09, there is a maximum in the 
| Vs 


horizontal displacement wz. In view of these 
results, we may say that the RayLEIGH wave 
appears at about #=3/. 


15 

(d) Against the results for the compres- 
sional origin, there are S phases in any of 
our figure for the distortional line source. 

(e) Putting ¢, 2>0o in (5.6)-(5.9), we get 
the following expressions for the permanent 
displacements at points far away from the 
origin 


— > 


ty af 5 (0). (6.2) 
Comparing (6.2) with (4.2), it seems rather 
strange that we have the permanent vertical (or 
horizontal) displacement larger than the hori- 
zontal (or vertical) one for the distortional (or 
compressional) wave origin. This question, 
however, may be settled if we consider that 
at v-oo, the radius vector from the wave 
origin to the observing point is almost parallel 
with the free surface. 

(f) The asymptotic values in (6.2) are shown 
in the figures for z=5/f and #=10/ in Fig. 
9. From these figures, we see that around 
these asymptotic values, the RAYLEIGH waves 
take their typical forms shown in Fig. 5, 


with # and v interchanged. 
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Propagation of Tremors over the Surface of an Elastic Solid. 
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and 
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Abstract 


Wave phenomena, which are seen over the surface of a semi-infinite elastic solid 
subjected to an impulsive line force of SV and SH type, have been studied. The results 


obtained for the SV type origin are almost the same as for the dilatational origin. 


In the 


case of SH type origin, there is neither P wave nor RAYLEIGH wave in seismograms to 


be obtained. 


§1. In a previous paper (H. Taxeucar and 
N. Kopayasut, 1954), under the same title as 
above, wave phenomena, which are seen along 
the surface of a semi-infinite elastic solid sub- 
jected to an impulsive force normal to the 
surface, were studied. In that paper, we got 
two formulus by which the surface displace- 
ments at any point and at any time can be 
calculated. By using these formulus, we could 
calculate the surface displacements in the 
neighbourhood of the wave origin and thus 
could make clear the way of appearence of 
RAYLEIGH waves which are the main parts of 
the displacements at points far away from 
the wave origin. In the present paper, we 
shall study the similar problem for a semi- 
infinite elastic solid subjected to an external 
_tangential force. The tangential force is 
assumed to be acting at («=0, y=0) in the 
rectangular coordinates shown in Fig. 1. There 
are two types of tangential forces acting at 
(c=0, y=0). In the so-called SV type of tan- 
gential wave origin, the force is acting parallel 
to the z axis. On the other hand, in the so- 
called SH type of tangential wave origin, the 
force is acting parallel to the zaxis. In 1904, 
H. Lamp studied the case of tagential wave 
origin of SV type and obtained the formal 
expressions for the displacements but he did 
not work out the results in detail. On the 


other hand, an experimental study for the 
tangential wave origin has recently been made 
by Nortawoop and Anperson (1953). In their 
paper, they said that when the force is of SV 


Fig. 1. 


type, the result is not so different from that 
for the normal wave origin. When the force 
is of SH type, however, the result obtained 
is entirely different from those for the above 
two cases. For example, there is neither P 
wave nor RAYLEIGH wave in the seismograms 
for the SH type wave origin. To study why 
these things occur will be one of the main 
objects of the present paper. 

§2. It was shown by H. Lamp (1904, his 
equation 55) that the horizontal and_ vertical 
displacement 2) and vp at (x, =O.ne)) dilemro 
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a SV type tangential force @ e’”’ acting at 
(e=0, y=0) are given by 


Q \ RV ee —k? einirén) ge | 


i —— Oe = Fe) 
Vog=t Q. ie £222 2 —h? V BR’) 
0 Tg esth exes FE) 
x elmirey de | 
where 
P(E) =(2&*—k?)? —4@* VP-P V2—k. 2.) 


In (2.1), the notations are all the same as in 
Lamb’s paper. In order to satisfy the radia- 
tion condition at 2 0 (A. SomMERFELD, 1949, 
p. 189), we must interpret the integrals in 
(2.1) as the limits of those along the path Z; 
in the complex € plane shown in Fig. 2 when 
semicircles around six singular points §=-h, 
-+-k and +« are made vanishingly small. We 
must also assume that //&—h? and VV &—R? 
are positive real on the right end of Z;. In 
order to calculate 2) and vp in (2.1), we shall 
consider 


uy=— 
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DoS Q Bel@t-K2) __ 7@) 
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ff €)=Q22—-P YP HEY PVP, 


(toe K2K2*—k?—2) e2—h? VV ne? — Rk?) 
F(t) 

where P denotes the principal value of the 

integral concerned. wy) and wv» in (2.3) are the 

formal expressions for the horizontal and ver- 

tical displacements due to the tangential force 


Qe, 
SO Lutting 
&=pab, ck=pc 


in (2.3) and generalizing the law of time varia- 
tion of the tangential force at the origin, we 
get 


h=pa, k=pb, (Sel) 


MT LUyp=P | U(A)QOt—adx) do, 
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Fig. 2 
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(2.2) 


x ei(ntt+ex) Z 


L3 in (2.2) is a large semi-circle with the center 


at 0. As there is no singular point inside 
(Z,+Z3), the integrals in (2.2) are equal 
to zero. Making the differences of the 


corresponding expressions in (2.1) and (2.2), 
we get 


“By In i k? 


i(pi-Ex) ‘ 
Lite 5 


et (nt- Ex) dé, 


(2a) 


MT pLVy = mMnBO(t—cx)+ \ V(A)Q(t—abx) dé 
1 
(ee) 


for the tangential force Q(Z) at (w=0, y=0). 
U(@) and V(@) in (3.2) are given respectively 
by 


4m? PV 0 —1 (m?—6") 


(262 —m?)* + 1664(@2—1)(m?—6?) 
Ue) = fOr el =< 0 ==777, | 
nal mV Pm an 
(26°—m?)?—40°/ 8 —1 V8 — me 
for m<6, 
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2020 — mV O—-1_ Vm? 


SONG: Cae 2. 92 
V6) = (26?—m?)' + 1664(62—1)(m?— 6?) 
for 1<é<m, 
0 for m<é, (ee 
where 
b Ve 1 
m= — = ) a= == p , 
Ws iene. y A+2n 
1 il 
Vs / Lt Vr Si 


Equation (3.2) is valid for any a(>0) and ¢. 
We shall now calculate zz) and vy in the case 

when 
Q(¢))=const Y, say for 
=() 


—T<fpxt 
(30) 


Puttin —ceeAnGl 7p—c2y) il (sez), we. have, 
after a few calculation 


otherwise. 
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—** 141) 

MT LUj=QVP 
Vie; 
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Vine ; 
: (FL) 
+@\ ; VO) do, (3.6) 
WEG 
—”"(¢!—1), 
xv 
where 6 is the Dirac’s function. This is the 


final formula by which we can calculate the 
“values of wz and wp at any 2(>0) and z 


§ 4. In our numerical examples, we assumed 


Vp 
Ve 1 3 


and calculated at first the values of U(#) and 
V(@) in (3.6). They are shown in Fig. 3 and 
4. Next, the values of les are U 
are calculated at = : =(). 1, 2 eine , Ine 


(4.1) 


b 
ely ff = 
a 


uM and Jo 


results are shown in Fig. 5 and 6. The results 
to be noted are as follows: 

(a) The horizontal displacement is outward 
almost everywhere. There is, however, a very 
narrow region of inward displacement. This 
region is propagated with the velocity of the 
RayLEIGH wave Vx=0.92Vs. The origin time 


U(e) | a 


Fig. 3 


Fig. 4 


of this wave is at ¢=—r. This is the time 
at which the tangential force begins to act on 
the otherwise free surface. 

(b) The fronts of both the horizontal and 
vertical displacements are propagated with the 
velocity of P wave Vp. The origin times of 
these waves are also at /=—tT. 

(c) In the horizontal displacements after 
¢=t, there appears a region of very large out- 
ward displacement. This region is propagated 
with the velocity Vr. The origin time of this 
wave is at f=r. This isthe time at which 
the tangential force is removed. 
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(d) The vertical displacement is upward in 
the domain between the wave front and the 
point (or rather line) at which the horizontal 
displacement is a large inward one, and it is 
downward in the domain between the points 
at which the horizontal displacements are 
large inward and outward ones respectively. 
Thus the domain of downward motion is pro- 
pagated with the velocity Vr. 

(e) The apparent wave length of the gener- 
ated RayLEIGH wave is about 2Vpr. This 
length is approximately equal to the extent of 
the downward domain at f=r at which the 
tangential force is removed. 

(f) The time variations of the displacements 
at points far away from the wave origin can 
approximately be obtained by tracing the a 
and wv» curves for ¢ =5 in Fig. 5 and 6 from the 


side of larger i to that of smaller a § 


Bigs Ve 


(g) The vertical displacement at 2=0 is 
zero all through the time. This will be under- 
stood from the consideration of symmetry. 
The horizontal displacement there is a sudden 
outward motion followed by a gradual decrease 
to the zero line. It is shown in Fig. 7. 

(hk) There is no S phase in any of our 


21 
figures. This result, together with some of 
the above results, is very similar to that 
obtained for the compressional wave origin 
already studied. 


§5. We shall now study the similar problem 
for an impulsive tangential origin of SH type. 
Refering to the rectangular coordinates in Fig. 
1, we have the following solutions of the 


equations of motion for a homogeneous elastic 
body. 


u=v=divi=0, 


w= Ae rttéx)-By | P?=&—PR?, 
Dd / 0 
b= Pa f 
Vg p 7D 
Tyz=T yy=0, 
en Ow = —pABei itt )-By 5 
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wy=0)=—+Ty. Y=0), (6.1) 
“B 


where A is a constant and the real part of 8 
is positive. Thus, for a tangential line source 
at «=y=0, we have 


wax, y=0, t)=Wo, say 


= Q aly d. 
Qu e Sey (8 o 
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Proceeding as in section 2, we can rewrite 
(G.2)eas 


is, a) oN elke - at etsr d | 
Biri Oe | shasta ese 
ee vay eo oes (5.3) 
T/A k p 


which, for the general type of time variation 
Q(t), is transformed into 


1 = O(t—bz8) 
o= 10, bee 
ve ft | 17 o°—1 ( ) 
where 
pas Fan BS een = POU (5.5) 
Vs 


In the case of an impulsive origin 


Q(é)=0 for t>d, 


this section. 
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STONELEY Waves Generated by Explosions. 


By 
Akira Kusorrera 


Geophysical Institute, Faculty of Science, Kyoto University, Kyoto. 


Abstract 


Using three components seismometers, field measurements of under ground motions 


within bore-holes were carried out. 


In this experiment, boundary waves which could be considered as STONELEY waves 


were obtained. 
9m depth from the ground surface. 


The boundary surface is between layers of clay and sand and is of about 


The characteristics of this wave are as follows: 


Tne wave velocity is 430 m/sec. 
The period of motion is about 0.06 sec. 


The particle trajectories are elliptical and retrograde. 
The amplitude of the wave is large at the boundary surface and is decreased rapidly 


on leaving this surface. 
These results are 
on the STONELEY waves. 


§ 1. Introduction 

Putting seismometers at various depths 
within bore-holes, and varying the depth of 
shot-hole, explosive waves were recorded by 
vertical and horizontal seismometers. 

The most interesting result obtained in these 
experiments was that the StoneLEY wave could 
be observed by explosions, at the boundary 
surface between two different layers. 

In 1924, R. Sronetey found that it was 
theoretically possible for RAYLEIGH-type waves 
to be transmitted along the boundary surface 
of two elastic solids. Later on, K. Srezawa 
and K. Kanat (1938 and 1939) studied mathe- 
matically the way of generation of the boun- 
dary waves. K. Srmzawa and K. Kanai (1939) 
and J. G. Scuourer (1947) calculated the pos- 
sible existence range of the STonELEY wave. 

Many theoretical studies on the STonELEY 
wave have been carried out since then. How- 
ever, there have been no reports on the direct 
observation of this wave or the opportunities 
of doing so. (Y. Sato 1953). 

In the present experiment, the seismograms 
recorded a phase which seems to correspond 
to the boundary wave. This phase could clear- 


in good agreement with those predicted by the classical 


theory 


ly be recorded on the trace of seismograms 
only when the seismometers were placed at 
or near the boundary surface of the two dif- 
ferent layers. 

In view of many other results obtained it 
seems that this wave corresponds to the 
STONELEY wave. 

§ 2. Experimental Procedure. 

Instruments 

The instruments employed in these experi- 
ments are those used in usual prospecting 
operations. The seismometers are of the 
E.T.L. type (Electro-Technical Laboratories’ 
Seismic Detector, Model EVS-3D) with three 
components (one vertical and two horizontal), 
the natural frequency being 10 cycle/sec., and 
the damping being critical. The galvano- 
meters have a natural frequency of 30 cycle/ 
sec., and are critically damped. Each galvano- 
meter is directly connected to a seismometer, 
with no amplifiers but sometimes with at- 
tenuators. 

Figure 1 shows a response curve of the 
entire system, including the seismometer and 
the galvanometer. The magnification of this 
system is about 2500 at 20 cycle/sec. 
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Magnification 


Frequency cycle/sec. ; 
Fig. 1. Frequency-magnification curve for Fig. 2. Geological profile of 
over all system including seismometer, the test field. 
galvonometer and occillograph. 
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Fig. 3. Time-distance curves by the refraction test. (after T. KAMETANI) 
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Location a 

The observations were carried out by the 
“ Seismic Exploration Group of Japan” in Feb. 
1954 and March 1955, on the ground of the 
agricultural field of Tokyo University, Tana- 
shi, Tokyo. 

The geological profile here is shown in figure 
2. The upper layer is a clay, called “ Kanto 
Loam,” under which is a sandy layer beginn- 
ing 9m below the surface. The StroneiEy 
wave was observed at the boundary surface 
between the clay and sand layers. 

A seismic prospecting test by the refraction 
method was carried out in this field by our 
group, and the result was reported by T. 
Kametani (1954). The time-distance curves in 
this report are copied in figure 3. The longi- 
tudinal wave velocity is gradually increased 
with the depth. This velocity distribution 
may be represented by the following formula: 


(1) 
where z is the depth in meter. Using the 


time-distance curves, the constants 7 and a 
are determined as follows: 


Vo=1411.5=-3.7 m/sec. 


v(Z)=Vo+az 


(= BVA seb ee 


Shooting procedure 

Two types of tests were made as follows: 
1) The ground motion resulting from dyna- 
mite explosions in a bore-hole of 15m depth 
was recorded by seismometers placed at bore- 
hole bottoms of varying depths. The depths 
of the seismometer holes were 2m, 7.5m, 
9.0m and 10m in the first test, and 1.5m, 


5.0m, 8.5m and 11.5m in the second test, 


the shot-receiver distances being 60 m and 45 m 
in the first and second tests, respectively. 
(Figure 4) 

2) In order to determine the propagation 
velocity of this wave, three component seismo- 
meters (vertical, horizontal radial and_ hori- 
zontal transversal components) were placed at 
three shot-receiver distances, namely 4=20m, 
4=60m and 4=100m along the boundary 
surface of the clay and sand layers, which is 
at about 9m depth. An additional instrument 
was placed at the up-hole 100m away. A 


Experiment 1) 


= 


iSm- 


Om- 


Sm- 


10m- 


15m-~- 


Fig. 4. 


Maps of disposition of seismometers 
and shot-holes in the experiment. 


dynamite charge was fired at the depth of 
15m under the ground surface. 

By using an attenuator, the sensitivity of 
each seismometer was kept constant in each 
experiment. The sensitivity of seismometers 
were changed, however, according to the 
distance as follows; 

Shot-Receiver Sensitivity Ratio 

Distance. Vertical Horizontal 
100 m 1 1 
60 m 1/3 il 
20m 1/10 1/10 


§ 3, Experimental Results. 

The seismograms obtained in experiments | 
and II are shown in figures 5, 6 and 7. 

As is seen in figures 5 and 6, there is a re- 
markable later phase in traces A-1, A-2, B-1 


and B-2, but no such phase in trace C, D, E-1 


and E-2. 

Our conclusion is that this clearly distinguish- 
able later phase is due to a boundary wave, 
on the ground that it is found near the boun- 
dary surface between the clay and sand lay- 
ers and its amplitude is decreased rapidly on 
leaving the surface. 

Is this wave progressive or not ? As shown 
in figure 7, we found in experiment II that 
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Seismograms obtained by the Ist test in 1) experiment. 


Shot-Receiver Sensitivity 
Distance Ratio Notation 
60 m 1 A-1 
Yi ” B-1 
YA Y E-1 
” " (e 
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Pigewos 
Trace Types of Disposed 
Number Seismometer Depth 
i Vertical 100/s. 10m 
2 y 9m 
3 y " PRR 
4 ” ” 73m 
Fig. 6. 
Trace Types of Disposed 
Number Seismometer Depth 
il Vertical 10c0/s. 14m 
2) i y 5m 
3 A Y 83 m 
4 " y 113m 


(Shot-hole depth; 15 m) 


this wave was recorded by the vertical and 
horizontal radial components of the seismo- 
meters at every distance on or near the bounda- 
ry surface, but not recorded by the up-hole 
seismometer. Thus, we may safely conclude 
that this wave propagated along the boundary 
surface. 

In the following paragraphs, some of the 
characteristics of these boundary waves will 
be studied. 

1) Propagation velocity 

The time-distance curves for the initial 
motions and boundary waves obtained from 
figure 7 are shown in figure 8. From these 
time-distance curves, we see that the propaga- 
tion velocities of the boundary wave and the 


Seismograms obtained by the 2nd test in 1) experiment. 


Shot-Receiver Sensitivity 
Distance Ratio Notation 
45m 1 E-2 
74 Va D 
w a iB=2 
" ” A-2 


initial motion are 430 m/sec. and 1800 m/sec., 
respectively. The propagation velocity of 
1800 m/sec. corresponds to the longitudinal 
wave by the following reason: 

By equation (1), the propagation velocities 
of the longitudinal wave at 9m (the boundary 
surface) and 15m depth (the depth of explo- 
sion point) are estimated to be 


v(2=9)=1728 m/sec. 
v(2= 15)=1920 m/sec. 
The velocity of 1800 m/sec. lies between 
these two values. 
2) Trajectories of Particle Motions. 
Figure 9 shows the particle motion trajector- 
ies obtained by the vertical and horizontal 


ie | ai 


‘STONELEY Waves Generated by Explosions. 


Fig. 7. Seismogram obtained by 2) experiment. 


Trace Types of Disposed 

Number Seismometer Depth 
1 Vertical 10c/s. 9m 

2 7 " . 9m 

3 ” ” 8m 

4 y Ua 0m 

15 Y 270 /s. Om 

6 Horizontal 1000/s. 9m 

if 1” y 9m 

8 ty y 84m 

9 1 y Om 
10 Y 2700/s. Om 


(Shot-hole depth; 15 m) 


radial seismometers at a point of 8.5m depth 
from the ground surface and 100m distance 
from the firing point. Three successive cycles 
are shown in the figure. The motion is re- 
trogade in the first two cycles and is roughly 
elliptical, but not so in the third cycle; It 
seems that the first two cycles correspond to 


the SvroneLEY wave but that the third one re- 


— 
ea 


Shot-Receiver Sensitivity 
Distance Ratio 
20m 1/10 
60 m aly es} 

-100 m it 
100 m 1 
100 m ; — 

20m 1/10 
60 m il 
100 m il 
100 m 1 
100 m — 


presents some other type of wave. 

At a point of 9m depth from the ground 
surface and 60m from the firing point two 
cycles can also be seen in trajectories. 

3) Pertod and Dispersion 

The period of this wave is about 0.06 sec., 
but. the time intervals between successive 
peaks and troughs become longer as time 
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Fig. 8. Time-distance curves for under 
ground test. 


goes on. This time increase may be due to 
the abnormal dispersion of the wave. The 
_STONELEY Wave is not dispersive, if the elastic 


Up Up 


constants and densities of the two layers are 
constants. Thus the wave dispersion obtained 
above may be explained by the complexities 
of the structure in this region. 

4) Amplitude. 

The amplitude distribution in a vertical 
column is shown in figure 10. On the other 
hand, the amplitude distribution versus shot- 
receiver distance along the boundary surface 
is shown in figurell. It is rather difficult to 
determine the attenuation formula of the 
STONELEY wave from this figure, because at 
4=20m, the SronELEY wave is not clearly 
separated from other waves. 

§ 4. Conclusions and Remarks 

In this section, the existence conditions and 
the excitation conditions of the StronELEY wave 
derived theoretically will be discussed. 

K. Sezawa and K. Kanai (1938) and F. G. 
ScHOLTE (1947) determined the range of the 
density ratio (o/o’) and the elastic constant 
ratio (#/#2’) of two separated layers within 
which SToNnELEY waves can be existed. Some 
of the results obtained are shown in figure 12. 
From figure 12, we see that the range of 


Up 


Fig. 9. Particle trajectories of succesive STONELEY wave recorded at 8.5m deep and 
100m shot-receiver distance. 
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Le oes ; 3 
ies A within which the StenELEY waves can 


be existed is quite narrow. When ( ai e~ 1) : 
fe] 0 


Ground 


Boundary 
Surface 


Surface 
1 


1 
1 
! 
\ 
3 ‘hd 
£ wer) 
ra 1 
£ ! 
4 i 
i 
t 
| 
I 
I 
0m 5m 10m 
Depth 
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Fig. 10. Effect of depth upon the vertical com- 


ponent of the amplitude of STONELEY waves. 


Aa 
20m 60m 100m 
Fig. 11. Effect of shot-receiver distance upon 


the vertical component of the amlitude of 
STONELEY waves. 


i.e., when the velocity difference of two 
separated layers is quite small, the StongLEy 
wave can be existed. 

In the present experiment, there was hardly 
any velocity difference between the two 
separated layers as was shown in figure 3, 
although the layers are clearly distinguishable 
in the bore-hole tests. 

Thus it may safely be concluded that in 
this field, the existence conditions for the 
STONELEY wave are satisfied. 

Next we shall consider the generation condi- 
tion of the StoneLEy wave: K. Sezawa and 
K. Kanar (1938 and 1939) showed mathemati- 
cally that dilatational wave origins or dis- 
tortional origins of SV type can generate the 


Ble’ 
3 
2 
1 
1 2 3 p/p! 
(Aje=2' |p =1) 
Pip 
3 
2 
1 
1 2 3 ple’ 
AQ=)'= co ) 
Fig. 12. The range of existence of STONELEY 
waves. (after K. SEZAWA & K. KANAI 


and F. G. SCHOLTE.) 
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boundary wave. In the present experiment, 
the origin is of dilatational type, because, in 
general, dynamite explosions initially generate 
spherical dilatational waves. 

The propagation velocity of the observed 
STONELEY Wave was measured as 430 m/sec. 
The propagation velocity of the SToNELEY 
wave is shown to be nearly equal to that of 
the transversal wave in the same solid. In 
the present experiments, the transversal wave 
velocity was not obtained, but assuming Pots- 
son’s ratio of about 0.45, we get the tranver- 
sal wave velocity of 430 m/sec. This value of 
0.45 is not so unnatural for clay or sand 
fields. 

Unfortunately, the attenuation formula for 
the StonELEY wave could not be determined 
experimentaly by the reasons given above. 
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